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Abstract 

Surface superconductivity is studied within Ginzburg-Landau theory for two 

classes of models for the order parameter of UPta. The hrst class assumes 

two independent one-dimensional order parameters {AB models), while the 

second assumes a single two-dimensional order parameter [E models). i?c3 

is calculated for all cases where the surface normal and magnetic held lie 

along high symmetry directions. Assuming specular reflection, it is found 

that except when H || c, the ratio Hc^/Hc2 is either unity or equals its 's- 

wave' value 1.695, although the precise H^z vs. T curve predicted by the AB 

and E models differs for various geometries. The results are compared with 

recent experiments, and predictions are made for future experiments. 
PACS numbers: 74.60.Ec, 74.70.Tx, 74.20.De, 74.60.-w 
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I. INTRODUCTION 



UPts is generally believed to be a reduced symmetry superconductor, i.e., one that 
breaks rotational and/or time reversal symmetries of the normal state, in addition to the 
U(l) gauge symmetry The strongest basis for this belief is the phase diagram which 

displays multiple superconducting phases just as liquid ^He has multiple superfiuid 

phases. In this paper we shall study surface superconductivity in UPta, and see if competing 
models make differing predictions which might be experimentally tested. 

We briefly review the existing models for the order parameter in UPts. There are two 
main classes. The first class, which we shall refer to as 'i?' models entails a single 

two-component order parameter, which may belong to any of the two-dimensional represen- 



tations {Eig, E2g, Eiu, E2u) of Dg^, the point group of normal UPts [|T4[. The second class 
of models posits two independent order parameters belonging to different representations 
T^ , [T5| -[T8| . In one subclass ||19|, which can be analyzed in great detail, comprising what we 



shall call ^AB^ models, both order parameters have the same parity; one transforms as an 
A representation and the other as a I? representation of D'qi^. Thus, if the common parity is 
even 



20| , the pair of order parameters transforms in one of four possible ways: (^ig, -Big), 



(y4ig, B2g), {A2g, Big), {A2g, B2g) ■ THs division into two classes is useful because in the 
simplest form of Ginzburg-Landau (GL) theory in which only terms that are formally of 
order (1 — T/Tc)^ are kept, all E models have the same formal GL free energy, and all AB 
models have the same GL free energy. This leads one to sometimes refer to them in the 
singular as the E model or the AB model. 

The E model invokes a symmetry breaking field (SBF), usually taken to be a weak 



antiferromagnetism that sets in at 5 K [|21| , |22| , in order to split the superconducting transition 
in zero field into two transitions as observed. The AB model posits two nearby transition 
temperatures by fiat. The original version of the E model was shown to be incompatible 
p!3| , p!5[| with experiments in that it fails to yield a tetracritical point in the H-T plane 



for H II c [E3|. ( Whether there is or is not a true tetracritical point in the experimentally 
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observed phase diagram is irrelevant. At the very least one has two phase boundaries Hc2{T) 
and Hci{T) that approach within 10-15 mK of one another, and the outermost line Hc2{T) 
has a sharp change in the sign of the curvature near this apparent tetracritical point. The 
theory can not explain this 'near collision' either.) This led us to propose the AB model 
which does not have this flaw. Our proposal in turn led Sauls to propose a refinement of the 
E model [||], which has three ingredients: (i) a specific order parameter (see below), (ii) 
a nearly cylindrically symmetric Fermi surface, and (iii) coupling of the gradient terms in 
the GL free energy to the SBF. A cylindrical Fermi surface causes the dangerous gradient 
coupling terms in the free energy ||23| to vanish within weak coupling BCS theory for an E2 
order parameter, and in conjunction with ingredient (iii) restores the tetracritical point for 
H II c. 

In this paper we shall use GL theory to calculate H^s for UPta within the original E 
model and the AB model. The GL free energy for the refined E model is mathematically 
isomorphic to, and can be regarded as a special case of, the AB model free energy . 



Since the key issue that determines H^s is the boundary condition on the order parameter, 
however, the refined E2U model of Sauls will have the same H^siT) curves as the unrefined 
E2U model for H ± c. Indeed, H^s curves may be similar for these two models even when 
H II c, except that one must beware of kinks in the refined model. We will limit ourselves 
to cases where the field H and the surface normal n lie in high crystal symmetry directions 
or planes. Further, we only consider ideal, or specularly reflecting surfaces. 

The motivation for this study is that we expect qualitative differences in the behavior of 
Hc3 between the various models for various geometries. Firstly, the boundary conditions are 
different amongst the E models, leading to differences in Hcs- Secondly, in the AB model 
the eigenvalue equations for ifcs decouple into separate equations for two components, and 
whenever the surface supports both components, the H^s curve is expected to mirror Hc2 
and show a kink. In the E model, on the other hand, even for H ± c, the gradient terms 
can couple the two components. It is thus possible for the kink in Hc3 to be smoothed out. 
We shall see that whether this happens or not is a question of dynamics, not symmetry. 



In fact, there turns out to be no smoothing with the gradient couphng values we use. A 
short quantitative estimate of this effect seems to be hard to get, however, so we present 
the full analysis which follows. Since an experimental measurement of ifcs has now been 
reported |^ for some of the geometries that we study, we can use our results to restrict the 
acceptable order parameters for UPts. We also hope that our work will spur a more detailed 
experimental study of Hcz in other geometries as well, as this will sharpen our understanding 
of the order parameter even further. 

We note here that as this paper was being written, we learned of a recent paper by 
Samokhin ||2^, also discussing ifcs in UPta within E models. Samokhin gives a microscopic 



foundation to the boundary conditions on the order parameter [^, while we take a purely 
phenomenological approach. Further, he focuses on the case where H || c, and the surface 
normal is arbitrarily oriented in the a-h plane. We have largely avoided detailed study of 
the H II c geometry because of the problems it presents in comparing to the observed bulk 
upper critical field. Thus Samokhin's work is nicely complementary to ours. 

The rest of the paper is organized as follows. In Sec. II, we recapitulate the bare 
essentials of the theory of surface superconductivity for fully symmetric superconductors 
with a complex scalar order parameter, paying special attention to the boundary conditions. 
We then extend these ideas to the AB and E models. The case of the E model is rather rich. 
Depending on the exact order parameter and field and surface orientation, the surface may 
or may not suppress superconductivity. We tabulate these cases and proceed to study them 
in Sees. Ill and IV. Certain technical aspects of the calculations are given in the Appendix. 



In Sec. V, we compare our results with the experiments of Keller et. al. ||2^, and see which 
order parameters are compatible. Our results are summarized in Table III, and the reader 
who is not interested in the details of the analysis should skip to that directly. We conclude 
with suggestions for future work. 

II. BASIC THEORY OF SURFACE SUPERCONDUCTIVITY 
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A. Fully symmetric superconductors 



We first recapitulate the theory for a superconductor with a complex scalar order pa- 
rameter ip, occupying the half space z > p8|-pO|l. Hcs is found by solving the linearized 



GL equation 

^(-^^V-|A)% = -a(T)^ (2.1) 

subject to the boundary condition of no current flow across the surface: j ■ z = 0. The 
notation in Eq. ( |2.1| ) is standard. 

If H II z, then then there is no surface superconductivity, and technically iJcs = -f^c2- 
The problem is more interesting when H ± z. For an interface with a vacuum or insulator, 
the boundary condition j ■ z = becomes 

D,^\,=o = 0. (2.2) 

where D is the gauge covariant derivative. Taking A = —Hzy, and writing ip as a plane 
wave in the x-y plane with wavevector {kx,ky) times a function f{z), Eqs. ( pTT|) and ( p.2|) 
can be rewritten as 



^2 f2nHY . ,2 ,2 



dz^ 



with {df/dz)z=o = 0. Here, zq = —ky^Q/27!-H , $o is the flux quantum, and ^(T) = 
—fi^/2'ma{T). The highest field eigenvalue is obtained with k^ = 0, and 

zo = goaT), H,, = -H,2 (2.4) 

with = 0.59010. [l/go = 1.6946.) Recall that Hc2 = hc/2ei^{T). 

Kittel's variational solution to the above problem is also worth reproducing. Taking 

/ = (2/7r(T2)i/4exp(-2V4(T2) (2.5) 

and minimizing the total energy with respect to zq and a gives 
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'2gK 

= (2/7r)i/V = 0.727^(T), (2.6) 

with gx = {I — 2/7r)^/^ = 0.603. {l/gx = 1-66.) Note that this variational solution has the 
property of zero integrated current in the y direction: 

/ jydz(x / {z-zo)fiz)dz = 0. (2.7) 
Jo Jo 

Finally, we note that if Eq. ( p.2| ) is replaced by the boundary condition ipi^) = 0, the 
highest H eigenvalue is obtained by letting — ^ oo in Eq. ( |2.3| ), and Hcs = Hc2- 



B. Reduced symmetry superconductors 

For both E and AB models, the order parameter can be written as a two-component 
complex vector 

This transforms as a single irreducible (E) representation of Dg;^ in the E model and as a 
reducible representation in the AB case. The quadratic part of the GL free energy density 
(which is all that matters for iJcs) can be written compactly as 

fSl = E MT)\vr\' + E ^:;{D.vrnD,vs). (2.9) 

r=l,2 i,j,r,s 

where i,j G {x, y, z}, and the k's are gradient coefficients suitably constrained by symmetry. 
It follows from Eq. (|2.9|) that the current is given by 

J^ = ^l^T.<JV:DJVs. (2.10) 

At first sight the condition j ■ n = cannot be written as a linear equation in rj. Just as 
in the case of ordinary superconductors, however, a linear equation can be obtained 



as was first done for ^He-A In contrast to ordinary superconductors, the roughness 
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of the interface is now very important. A rough surface is in general pair breaking. Since 
we wish to study the maximum possible enhancement of H^^, we will only consider ideal 
or specularly reflecting surfaces. Further, as mentioned before, we shall only consider cases 
where the surface normal n and the field H lie in high symmetry directions or planes. Under 
these conditions, the microscopic analyses P^ , |27|pT| , which we shall not repeat, show that 



the boundary condition on each component of rj reduces to |32 



r/^ = or (n-D)r/,. = 0. (2.11) 

Note that both components need not obey the same condition; one could have r/i = and 
(fi- D)?72 = 0. The general rule can be written as follows. Let us write the momentum-space 
gap function in the Balian-Werthamer notation, 

2 

Ao(k) = ^?7r^r(k), (even parity), (2.12) 

r=l 
2 

(i(k) = ^ r7r'?Ar(k), (odd parity), (2.13) 

r=l 

— * 

where ipr or ip^ are basis functions of the appropriate representation. Note that in the odd 
parity case, d and ifj^ are vectors in pseudo-spin space. The precise boundary conditions in 
Eq. ( p.ll|) are then as follows. Let us denote by Rk the wave vector obtained from k by 



reflection in the plane normal to fi. Then, in the even parity case, at the surface we have, 

r]r = Q if ^^(Rk) = -V^^(k), (2.14) 
(fi ■ T>)r]r = if V^^(Rk) = +^r(k). (2.15) 

The surface acts in effect as a momentum space mirror. If this mirror plane is a nodal plane 
for ipj., then tpy. must vanish at the surface. If the mirror plane is antinodal, then the normal 



derivative must vanish. In the odd parity case, we replace ipr by ipr- As noted in Ref. |^ 
it is unlikely that in reality ipr will be exactly odd or even under reflection normal to any 
n, even those of high symmetry, as the surface will in general flip pseudospin. Keeping our 
goal of studying surface superconductivity under idealized or optimal conditions in mind, 
we shall ignore this complication. 



In Table I we give a few illustrative cases of boundary conditions obtained by applying 
the above rule to candidate order parameters for UPts. Since we will always take H ■ n = 0, 
it is always possible to choose A ± fi, and the boundary condition ( p.l5[ ) reduces to VnVr = 
r]'^ = 0. Several points about this table should be noted: (i) The directions x, y, z are fixed 
along the crystal symmetry directions with x = a, y = a*, and z = c. (ii) For odd parity, 
unit vectors in pseudospin space are denoted by a subscript 's'. (iii) For each representation 
we have generally given the simplest possible basis function. As noted in Ref. this does 



not yield the most general gap function allowed by symmetry. As a counter-example, we 



show two order parameters for [34|, distinguished by superscripts. The second of these 



is in fact the one advanced by Sauls [jl| and by Norman |3^. (iv) In each case, we list the 



pair of quantities that must vanish at the interface. Thus for the example, with n = x, 
we must have f]'i = f]2 = 0. (v) In some cases, the boundary conditions are mixed, i.e. of 
the type discussed in Refs. [ p6| , ^ . 

The reader can write down the correct boundary conditions for any other order parameter 
using our rules. In subsequent sections, we will solve the H^s problem for both AB and E- 
type models. In the case of the AB models, H^s is independently found for the A and B 
component. By a simple rescaling of the coordinates, the eigenvalue equation can be reduced 
to that of the isotropic case. The boundary condition is either 77^ = or 77^ = 0, yielding 
Hf.-i/Hf.2 equal to 1 or ^ respectively. Having found Hf.-i{T) in this manner for each 
component separately, the thermodynamic Hc^{T) is taken to be the larger of the two for 
each value of T. A similar procedure will turn out to work in many cases for the i?-models, 
although it is not obvious a priori that the equations for rji and 772 will decouple. 

III. i?c3 FOR AB MODELS 

In this section we will calculate ifcs for the AB model order parameters listed in Table 
I. Writing {rja.rjb) instead of (771,772), the free energy (to quadratic order in rf) for all AB 
models is given by 
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/gL = [ar|?7r-P + /tr|D_L?7r-P + /tr|D2?7rP ■ (3.1) 

r=a,b 

Here Ur = ao(T — Tr), and all k's are positive. We assume that rja and rjb belong to the A 
and B representations, respectively. 

A. n = z 

Since all directions of H _L n are equivalent, let us take H = HSt, A = —Hzy. Mini- 
mization of /ql leads to the GL equations, which decouple for rja and rjb: 



52 f d 2ie ^ 



92 

Vr - '^'r-g^Vr = -arVr- (3.2) 



dx'^ \ dy he 
The boundary conditions are dzTja = rjb = at z = 0. 

Equation ( p.2| ) can clearly be made isotropic by a rescaling of coordinates. We then get 

H^^{T) = im5H:,{T), (3.3) 
Hl,{T)=Hl,{T), (3.4) 

where the superscript labels the order parameter component that is nucleating. 

Taking the Hc2{T) curves to have a kink as seen experimentally, we obtain a very simple 
picture for Hcs- For the even parity case, if > T^, there is only one branch to Hc3{T), 
and it onsets at T^, i.e., the upper transition in zero field. If > Ta, there is again only 
one branch, but it onsets below the upper transition temperature, and it extrapolates to 
for H = 0. These possibilities are shown schematically in Figs. 1 and 2. For the odd parity 
cases, the roles of rja and rji, are reversed, i.e.. Fig. 1 applies if Ti, > Ta, and Fig. 2 applies if 

Ta > n. 

B. n = X or n = y 

Let us first consider the Aig © i?2g order parameter. If we take fi = x, H = Hy, the 
analysis is identical to that just given, and we recover Figs. 1 and 2 for Ta > T^ and > Ta 
respectively. If n = y, and H = i/x, say, then taking A = Hyz, the GL equation becomes 



/ 9^ d'^ \ ( d 2ie Y 

which can be isotropized as before, but the boundary conditions now are dzfja = dzTjb = 
at z = 0, so Hcz is given by 

H:^{T) = l.mbH:^{T), r = a,b. (3.6) 

The resuhing picture is shown in Fig. 3. Note that within the GL approximation, the H^s 
and Hc2 kinks occur at the same temperature, as can be proved by simple geometry. Further, 
in this case, we expect to see inner transition hues analogous to those for Hc2, corresponding 
to the onset of surface superconductivity in the other component. 

We get the same answers as above for the Aiu © Biu the boundary conditions 

are the same. For the Aiu © order parameter, the results for the fi = x and the n = y 
cases are interchanged. 



IV. Hc3 FOR E MODELS 

We turn at last to our main problem, namely that of Hcs for the E models. For the 
reasons mentioned in Sec. I, we will not consider explicitly the refinements made by Sauls 

i- 

The free energy for all E models is more conveniently represented by regarding (?7i,?72) 
as the components of a vector in the x, y plane. This is self-evident for the Ei cases, but as 
shown by Tokuyasu , the number and form of the invariants is the same for the E2 cases 
as well. Writing {rji,rj2) as {rj^^riy), we can write the quadratic part of /gl as 

/gL = (^AVy? + (^AVx? + KlD*T]*DiT]j + K2D*T]*DjT]j + KsD*T]*Djr]i + K^DlT]*DzT]i, (4.1) 

where, i,j G {x,y}, and we sum over repeated indices. Also, we have a± = ao(T — T±), 
with T± = Tco ± e/ao. 

We write the GL equations implied by Eq. ( [4.1] ) using the notation of Ref. We 
scale all lengths by I = {hc/2eHy/'^, a^T and e by k;,//^ with k^, = ki + (K23/2), and 
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GL 



(4.3) 



/^23 = ^2 + Ks. We also define pi = —iDi, u = K23/2k;,, v = {k,2 — k^)/2ki,, w = k^/ ki,. Then, 
the GL equation is 

\'ly / \ 'ly / 

with 

(1 + u)pI + (1 - u)pI + wpl + e u{p^,py] + ivHjH 

u{Px,Py} + ivHjH {1 - u)pI + {1 + u)pI + wpl - e_ 

Here, {px,Py} is an anticommutator, e = eP/nb, and = ao(T— Tco)/^/^^. Note in particular 
that i (xl/H. 

To avoid misunderstanding, we emphasize that in this paper we will always take the axes 
X, and z as fixed along the crystal axes a, a*(± a), and c. This convention differs from 
that adopted by many other workers who are guided by the near isotropy of Hc2 for basal 
plane fields p5|j38| -^. This effect is generally explained by assuming that the staggered 



magnetization M''" rotates so as to be orthogonal to H. It is then convenient to rotate the 
axes, and the pair [rj^^riy) so that H is always along x (or y), and M"!" is implicitly along y 
(or x). We shall not do this. Rather, we shall rotate H and M''", and use a fixed fc-space 
form for rjx and rjy. The coupling between and rj then changes as M"!" rotates, and for 
M''" along the x or y axes, in particular, it takes the form 

r-C|Mt|2(|ry,|2-|,^j2)^ Mt||x, 

FsBF = <^ (4.4) 
l-C|Mtp(|ryj2_|^^,|2)^ Mt||y, 

where C is a coupling constant. Thus, strictly speaking, with our convention, Eq. (|4.3| ) is 
correct as written only if H lies in the x-z or y-z planes, and the sign of e depends on 
that of C. The latter depends on presently unknown microscopic physics. It can be seen, 
however, that Eq. ( |4.3| ) has a symmetry with respect to a change in sign of u. By taking 
sgn(M) = sgn(C), we can explain the observed kink in H(.2i whatever sgn(C) may be. 

To calculate ifcs, we must add boundary conditions. Before discussing these, however, 
let us comment further on the parameters m, f , and w. Firstly, local stability of the uniform 



solution requires that [|10 
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w; > 0, l-v>0, l + v>2\u\. (4.5) 

It follows from these that \u\ < 1. Secondly, the H ± c case involves w and u, but only \u\ 
can be found from experimental Hc2 slopes. (The H || c case involves only v, but comparison 
with experiment is problematic as discussed in Section I.) The values of u so determined 
vary quite a bit. We shall take 1^1 = 0.46 in our numerical work in accord with the data of 
Ref. §. 



A. n = z, H _L z 



Because of basal plane isotropy, we expect the same Hc3 curve for all field orientations 

in this case. Let us take H || x, e > 0, and let us work in the gauge A = —Hzy. Then 

in Eq. (^^ ), = —idzi Px = —id^, Py = —idy + z. We can clearly take ?7 as a plane wave 

in the x-y plane and replace by k and "Py by Zq — z, where k and —Zq are wavevector 

components, Then, {pxiPy} 2k{z — zq), and we get a one- dimensional Hamiltonian, 

-wd^ + {1 ~ u){z - zof + {I + u)k'^ + e 2ku{z - zq) 

2ku{z - zo) -wdl + (1 + u){z - zof + (1 - u)k'^ - I 

(4.6) 

For the £'2^, E\ui and e'2u cases listed in Table I, the boundary conditions are d^rix = 
dzVy = at 2; = 0. If = for all T, we get two decoupled ifcs problems and a kink in H^s 
as in Fig. 3. The only nontrivial features arise from the possibility of having k ^ 0. Such 
/c-mixing is found to affect the Hc2 curves for u > 0.61 for basal plane fields Whether 



GL 



the lowest energy eigenvalue for the Hc3 problem ever corresponds to /c 7^ is a dynamical 
question, to which we have not been able to find a two-line answer. It is obvious that having 
k ^ raises the energy from the diagonal terms in Eq. ( [4.6| ). The compensation of this 
rise from the off-diagonal terms is likely to be largest when the x and y solutions are nearly 
degenerate, i.e. for fields near the kink found with A; = 0. If the energy were lowered by 
having k ^ 0, the crossing of the x and y solutions would turn into an anticrossing, and the 
kink would be smoothed out. Our goal is to explore this possibility. 
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It is clear from Eq. ( [4.6|) , that a k ^ solution is more favored for large values of u. 
Let us mention two obvious analytic approaches to this problem. Firstly, a lower bound on 
the minimum u required can be found by constructing a Schwarz inequality for the matrix 



different widths and ay for rjx and riy. The energy is then a function of five variational 
parameters, a^, Cy, zo, k, and a mixing angle, and the minimization is complicated enough 
that one might as well attempt a nonvariational numerical solution. This approach does 
reveal one fact, however. If u is small, ax and ay are likely to be similar, and Eq. ( p.7|) shows 
that the matrix element of {z — Zq) is likely to be small. Thus the effect of the off-diagonal 
terms in Eq. ([4.6|) is probably of O(m^). 

We have therefore solved the eigenvalue problem for if^s numerically. We do this by 
working in a basis of harmonic oscillator eigenfunctions with unequal length scales (w/(l =F 
u)Y^^ for r]x and r]y. Only even parity (about z = 0) functions are used in accord with the 
boundary conditions. The Hamiltonian ( |4.(j| ) is almost diagonal in this basis. The evaluation 
of the off-diagonal terms is described in the Appendix. We find the lowest eigenvalue E^in 
and associated eigenfunction Ir^min) of i^cL numerically for given k and 2:05 and minimize 
with respect to k and zq via the conjugate gradient method. The required gradients are 
efficiently found using the identity 



and a similar one for the zq derivative. These calculations are done for e ranging from 0.1 
to 10. This range is large enough to encompass the kink region comfortably, as all terms in 



We find that for u = 0.46, k is numerically equal to zero for all e values considered. In 
fact, fc-mixing occurs only for u > 0.59, but the resulting anticrossing is very narrow, and 
substantial curvature in Hcz{T) is only visibly evident for u > 0.8. Thus, for experimentally 
relevant values of u, the picture for Hcs is the same as in Fig. 3, i.e., the Hcz/Hc2 ratio is 



element of the off-diagonal term in Eq. ( [4.6[ ), but this bound is too small to be useful. 
Secondly, one might try and construct variational solutions analogous to Eq. (|2.5|), with 




(4.7) 



Eq. ( [4.6| ) are of order unity. 
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always ideal and there is a kink in Hc3 Once again, we expect to see inner Hcs lines 
which we have not explicitly calculated. 

The Hc3 problem is trivial for the Eig and the e!^ cases listed in Table I. The boundary 
conditions now are rix = rjy = at z = 0, and the surface cannot support superconductivity, 
i.e., = Hc2- 



B. n _L z, H II z X fi 

In this geometry, as can be seen from Table I, except for the e!^J order parameter, the 
boundary conditions call for only one of the rji to vanish at the surface. The other component 
is not suppressed by the surface. The form of Hc3{T) is different depending on whether the 
surface suppressed component has the higher or lower transition temperature. This in turn 
depends on the sign of the coupling of the symmetry breaking field (SBF) to the order 
parameter. 

To understand this point, let us first consider the Eig and Ei^ cases from Table I. The 
two high-symmetry geometries are shown in Fig. 4. For cases (a) and (b), we have M''" || x 
or M''" II y, and e = — C|M'''p or e = +C|M'l'p, respectively. Suppose C > 0. Then, for the 
case of Fig. 4(a) we have rj'y = and e < 0, i.e., rjy has the lower T^. For the case of Fig. 
4(b) we have = and e > 0, i.e., has the lower T^. In both geometries, the surface 
and the SBF act in opposition, i.e., the SBF is oriented so that the surface supported order 
parameter component has the lower Tc. In contrast, if C < 0, then the surface and the SBF 
act in concert, and the surface supported component has higher Tc. 

The sign of ( is unknown and fixed by microscopic physics. We can conclude, though, 
that whether the surface and SBF act in concert or opposition, for Ei order parameters, the 
same behavior will be seen for the two geometries in Fig. 4. This point can be understood 
as follows. An Ei order parameter behaves as a vector in the x-y plane, and so if we rotate 
the surface normal n, keeping H || z x n, the boundary conditions on rj rotate in the same 
way, and we get an isotropic i^cs- 
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For the E2g and E2J order parameters, we can see from Table I that the same component, 
rjx, is supported by the surface for both n = x and n = y. For any sgn(C), this component 
must have the lower for one orientation, and the higher for the other, leading to 
different answers for Hcs- In terms of symmetry, this can be understood by saying that as 
n is rotated in the basal plane, an E2 order parameter rotates at twice the rate of a vector. 
Denoting the angle between n and x by 0, the boundary conditions change with (f), reversing 
every 30°. Thus the SBF and the surface will act in concert for one geometry in Fig. 4, and 
opposition for the other. 

We study the cases of surface and SBF in concert and in opposition separately below. 
It is convenient to fix e > 0, so that riy always has the higher Tc, and to rotate n and H as 
needed. 



GL 



1. Surface and SBF in concert 

We take the geometry of Fig. 4(a), and rj'y = rj^ = at x = 0. We also take A = —Hxz, 

allowing us to put pz = {—idz + x), px = —idx, Py = —idy in Eq. ( [4.3| ). Further taking 77 

as a plane wave in the yz plane, we can replace Py by k, pz by (x — Xq), and {px,Py} by 

—2ikdx- (Obviously, k and xq are wavevector components.) The resulting one dimensional 

GL Hamiltonian is 

-{I + u)dl + {I — u)k'^ + w{x — xqY + e —2iukdx 

—2iukdx ^(1 ^ tJ-)dx + (1 + u)k'^ + w{x — Xq^ — e_ 

(4.8) 

If = 0, we get two decoupled problems, and because of the boundary conditions 
if^3 = 1.69if^2) — -^c2- This leads to the phase diagram of Fig. 1. 

Thus, as in subsection A, the only real question is whether k gives a lower energy. 
We have diagonalized i^cL numerically using a harmonic oscillator wavefunction basis with 
length scales [{l±u)/wY^^, and odd and even parities about x = 0, for r]x and r]y respectively. 
The calculation is done as in subsection A, and we again find that with u = —0.46, k = 

15 



for all values of e [H2|. The Hcz picture is that of Fig. 1. 



2. Surface and SBF in opposition 

We take the geometry of Fig. 4(b), and rjy = rj'^ = aX y = 0. Further, with A = Hyi, 
and a plane wave dependence in the xz plane for r/, we get a one-dimensional problem in 



the y direction. Numerical work again reveals that A; = for all e so the i^cs picture is 
as shown in Fig. 2. 

C. n _L z, H II z. 

The high-symmetry geometries are shown in Fig. 5. In this case the field can not orient 
the SBF. Let us first assume that M''" stays fixed parallel to a*, i.e., that the appearance 
of superconductivity at the surface does not reorient M'''. Now, by the argument of the 
previous subsection, we expect the same i^cs behavior in the two geometries for the and 
i?2M order parameters, and opposite behavior for the the Ei order parameters. 

The two types of behavior, i.e., surface and SBF in concert and opposition, can be seen 
by studying the Eig case with ^ < 0. For Fig. 5(a) and 5(b), we have rj^ = rj'y = Q aX 
a: = 0, and r^^ = ?7y = at y = respectively. Since H || z, there is no kink in Hc2 and 
none is expected in Hcz either. The ratio Hc3/Hc2 is still of interest and we have found this 
numerically. When rj^ = 0, [Fig. 5(a)], surface superconductivity is not supported in the 
higher Tc component, and we find, as expected, that Hcs = Hc2- When rj'^ = 0, [Fig. 5(b)], 
surface superconductivity is supported in the component that onsets first as T is lowered, 
and > Hc2- Our results with u = 0.46, and v = 0.1 are shown in Fig. 6. The ratio 
Hcz/ Hc2 varies from 1.71 for T Tc+ to 1.77 for T Tc+ — 2ATc. 

The assumption of fixed M''' parallel to a* may not be physically relevant, however. Since 
the observed magnitude of M''" is very small [^,0 , reorientation of M"'" near the surface must 



be considered. If || a near the surface, and || a* in the bulk, the superconducting 
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condensation energy could outweigh the magnetic anisotropy and domain wall energies. In 
this case, we expect a picture resembling Fig. 6 for all the E models for any fi _L z, H || z. 

V. COMPARISON WITH EXPERIMENT, PREDICTIONS, AND CONCLUSIONS 



As stated before, Keller et al. have reported surface superconductivity in a whisker of 
UPta. The long axis of the whisker is along the c-axis of the crystal, and it has six facets along 
the {100} crystal planes, i.e., the surface normals are n = x = a, and five others obtained 
by successive 60° rotations about z. Keller et al. study the ac resistivity and susceptibility 
for various field orientations, two of which are H || a* and H || c, corresponding to the 
geometries in Figs. 4(a) and 5(a) of our paper. The magnitude of the ratio of the critical 
field to the bulk Hc2 (1.19 and 1.35 for H || a*), plus the cusp-shaped deviations in the 
critical field for fields with a small (positive or negative) component along n (see Figs. 1 
and 2 of their paper), leave little doubt that the observation is indeed of ifc3- 

Let us begin by studying the case where H || a* = y. Recall that fi = x. Keller et al. 
report Hc3/Hc2 = 1-35 for T = 470 mK. Given that (i) Tc+ = 508 mK, and that the bare 
Tc splitting 2e/aQ, is about 35 mK in the annealed crystals, and (ii) that Hc3/Hc2 decreases 
to 1.15 and 1.05 for T = 450 and 415 mK, respectively, which are above and below the 
tetracritical temperature, T* = 430 mK, we tentatively conclude that the H^s curve is of 
the type shown in Fig. 1, i.e., there is surface superconductivity in only one component, 
and it onsets at Tc+. This picture is obtained within the Aig © i?2g and Aiu © Biu models 
if Ta > Tb. Within the Aiu © model the predicted picture is that of Fig. 3, whether 
Ta > Tf, or Tf, > Ta- This does not seem to fit the observations, which would appear to 
exclude the specific Aiu © B2u order parameter listed in Table I. Within the E models, the 
picture requires that rjx have the higher Tc when H || y, i.e. || x. This in turn requires 
C < for the Ei models, and C > for the E2 models. These facts are summarized in Table 
II. 

Let us now consider the H 11 c, n = x case. An H^s > Hc2 is seen at T = 0.932Tc ^ 473 
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mK. Since data are not given at any other T, the picture of Hc^{T) could resemble any of 
our Figs. 1, 3, or 6. For the AB models we expect the same qualitative behavior as when 
H II a*, since the direction of H in the plane normal to fi does not alter the boundary 
conditions. For any of the four E models (-Big, E2g, Eiu, E2u), on the other hand, a little 
work shows that if we assume that M''^ stays fixed along a*, then the sign of ( in Table II 
is such that for the present geometry the surface and SBF are in opposition. This would 
imply no surface superconductivity. As argued in Sec. IV C, however, a reorientation of 
is quite possible, in which case all four E models would lead to an Hcs curve as in Fig. 6. 

We thus see that the present ifcs data can only be used to exclude (and not very con- 
fidently at that) the Aiu ® order parameter in Table I. We can however predict the 
expected behavior for the other geometries based on the results in Sees. Ill and IV, given 
the restrictions in Table 11. These predictions are summarized in Table III, which is the main 
result of our paper. 

It can be seen that widely different results are obtained for the geometries n = c and 
n = a*. The case fi = c distinguishes between Ei and E2 models and between E and AB 
models, if one can confidently detect the presence or absence of a kink in ifc3- The case 
with n = a* and H || a distinguishes between the odd and even parity E models, and the 
AB models. Systematic experimental efforts to study these cases, would be, to state the 
obvious, extremely valuable. Even for the one surface that is presently available, fi = a, 
studying the variation of i^cs with field orientation in the a*-c plane would be useful and 
would impel more theoretical study. 
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APPENDIX: 



We describe here our numerical procedure for finding the matrix of the GL Hamiltonian 
(|4.3|), using the case of Sec. IV A for specificity. We denote the operator valued elements of 
the matrix ( |4.6|) by H^x, H^y, Hy^, and Hyy, and define 

H^^^ = -wdl + {It u)z\ z = x,y. (Al) 

We work in the eigenbasis \n,i) of H^p: 

Hl-^\n,i) = [w{lTu)f^ i2n + l)\n,i), i = x,y; n = 0,1,2, ... . (A2) 

These functions are normalized on the half-space z > 0. Depending on the boundary 
conditions, only even or odd n may be needed. 

The matrix representation of ifcL in this basis is clearly very simple, except for the 
following elements: 

Pnm = {n,x\z\m,x), (A3) 

Qnm = {n,x\Tn,y), (A4) 

Rnm = {n,x\z\m,y). (A5) 

We do not need to define {n,y\z\m,y) separately, as 

(n, y\z\m, y) = [(1 + u)/{l - u)Y'^Pnm. (A6) 

These matrix elements can be found efficiently using the generating function of the 
Hermite polynomials. We show how this is done for Rnm- Note first, that 

Rnm = ^ = / zexpl-^-^]Hnl^]H^\^]dz, (AT) 

^7r^x^y2^+^n\m\ V J J \^y J 

where ^x,y = [""^/(l -F and if„ is a Hermite polynomial. Denoting the integral in 

Eq. ( [A7| ) by Inm, and using the generating function, we get 
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(A8) 



2A 



^{B'^-4AC)/4A 



— =erfc f — = I + erfc' f — j= 
2^fA \2^fA) \2^fA, 



where A = + ^y^)/2, B = -2{s/^^ + t/Q, C = + t^, and erfc'(x) = deifc/dx, with 

, , , 2 / x"^ \ , . , 

erfc(.) = l-_ + + . (A9) 



The problem of finding Rnm is now reduced to expanding the right hand side of Eq. (|A8|) 
in powers of s and t, and comparing coefficients. This can be done as follows. Consider a 
polynomial A{s,t), of maximum order s^t^^: 

N M 
ri=0 m=0 

The array defines the polynomial. The array for the sum of two polynomials A{s, t) 
and B{s,t) with arrays anm and bnm is obviously 

Cnm = anm + bnm, < U < N, < 171 < M. (AH) 

The product A{s,t)B{s,t) has terms up to order s^^t^*^. If we are only interested in terms 
up to order s^t*^, however, we can write the corresponding array as 

n m 

Cnm = Y.Y. aijbn-i,m-j, < U < N, < TU < M. (A12) 

i=0 j=0 

Note that only array elements anm and bnm with n < N, m < M appear in this formula. 

The above procedure can be used to expand the right hand side of Eq. ([A8| ) . Suppose 
we want for n,m < N*, where A^* is a number of order 20, say. We write B and C 
as polynomials in s and t, and represent them by arrays as above. ( Most of the elements 
of these arrays are zero at this stage.) We then Taylor expand exp{B^ /AA), exp(— C), 
eTic{B / 2^/A) , and erfc'(i?/2-\/]4), and repeatedly use the array operations ( |A11| ) and ( |A12| ) 
with N = M = N*, to evaluate Eq. (|A8D . The Inm can be directly read off the resulting 
array. 
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FIGURES 

FIG. 1. Schematic behavior of Hcz{T) when only the higher Tc order-parameter component 
supports surface superconductivity. 



FIG. 2. Same as Fig. 1, when only the lower Tc component supports surface superconductivity. 

FIG. 3. Same as Fig. 1, when both components support surface superconductivity. In this case 
we expect inner H^z lines (dot-dashed) near the kink in 11^^. We have not explicitly calculated 
these lines. 



FIG. 4. High symmetry geometries when n _L z, H || z x n. a, a*, and c are the crystal axes. 



FIG. 5. Same as Fig. 4 for n _L z, H || z. 

FIG. 6. Behavior of Hc^{T) in the E models for n _L z, H || z, when the higher Tc component 
supports surface superconductivity. The dot-dashed line is the ratio Hc2,/Hc2- The curves shown 
are calculated for u = 0.46, v = 0.1. 
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TABLES 

TABLE I. Boundary conditions for various candidate order parameters for UPta. The last 

three cohimns hst the pair of quantities that must vanish at the surface for given surface normal 
n. Note that (i) a prime denotes a normal derivative, (ii) the axes x, y, and z are fixed along the 
crystal symmetry axes a, a*, and c, respectively. 









Quantities that vanish for 




Representation 


tpr or Ipr 


n = z 


n = X 


n = y 


Elg 






{Vl,V2) 


{v'i,m) 


E-2g 


i^x ~ "^kxky) 








Elu 


{kxj ky^Zg 


{V'1,V'2) 


{m,V2) 


iv'1,112) 




i^kxXg kyl/g^ kxUs ~t~ ky'^s) 


{V'1,V'2) 


Mixed 


Mixed 




{kx — ky,2kxky)kzZs 


{vi,m) 




{Vl,V2) 


Alg e B2g 


(1, kxkz — Skxkykz) 






iVl,V2) 




{kz, ^kxky — ky)zs 


{Vl,V2) 


{V'1,V'2) 


iVl,V2) 




{kz,3kxky — k^)zs 




{v'i,m) 


Wl^V2) 



TABLE IL Constraints on the transition temperatures, or the sign of the coupling to the 

symmetry breaking field, for the order parameters listed in Table I, obtained by requiring agreement 
with the Hcz data of Keller et al. for H || a*, n || a. 

Order parameter Constraint 

Alg © B2g Ta > n 

Alu © B2U No agreement 

Am © Bm Ta > n 

Elg or Elu C < 

E2,, or 4^,^ C > 
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TABLE III. Expected iJcS vs. T behavior for the principal geometries. The order parameters 
for each case are as in Table I, and the signs of {Ta — Tf,) and ( are as in Table II. A blank denotes 
a case that we have not studied. 'No SSC stands for no surface superconductivity. 



Order parameter 


fi = c 
H±c 


fi = a 

H II a* 


fi = 

H II 


a* 
a 


fi = 
H II 


a 
c 


fi = a* 
H II c 


Elg 

E2g 

Elu 

E2u 
^2u 


No SSC 
Fig. 3 
Fig. 3 
Fig. 3 

No SSC 


Fig. 1 
Fig. 1 
Fig. 1 

Fig. 1 


Fig. 
Fig. 
Fig. 

Fig. 


1 
2 

, 1 

2 


Fig. 
Fig. 
Fig. 

Fig. 


6" 

ga 
ga 

ga 


Fig. 6 
Fig. 6" 
Fig. 6 

Fig. 6" 



Aig ® B2g Fig. 1 Fig. 1 Fig. 3 Fig. 1 Fig. 3 

Aiu ® Fig. 2^ Fig. 3 Fig. F Fig. 3 Fig. F 

Aiu ^ Bu, Fig. 2 Fig. 1 Fig. 3 Fig. 1 Fig. 3 

a. Assuming that the surface reorients Otherwise no surface superconductivity is 
expected. 

b. Assuming > 7],. Otherwise Fig. 1 apphes. 

c. Assuming Ta > T^. Otherwise Fig. 2 apphes. 
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